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Abstract. Be/X-ray binaries are systems formed by a massive Be star and a mag- 
netized neutron star, usually in an eccentric orbit. The Be star has strong equatorial 
winds occasionally forming a circumstellar disk. When the neutron star intersects 
the disk the accretion rate dramatically increases and a transient accretion disk 
can be formed around the compact object. This disk can last longer than a single 
orbit in the case of major outbursts. If the disk rotates faster than the neutron 
star, the Cheng-Ruderman mechanism can produce a current of relativistic protons 
that would impact onto the disk surface, producing gamma-rays from neutral pion 
decays and initiating electromagnetic cascades inside the disk. In this paper we 
present calculations of the evolution of the disk parameters during both major and 
minor X-ray events, and we discuss the generation of gamma-ray emission at different 
energies within a variety of models that include both screened and unscreened disks. 
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1. Introduction 



A significant number of variable gamma-ray sources of likely galactic 
origin has been found by satellite-borne experiments such as COMP- 
TEL and EGRET (e.g. Torres et al. 2001, Zhang et al. 2002, Nolan 
et al. 2003). These sources vary their flux on timescales from weeks 
to months, similar to the timescales exhibited by X-ray transients at 
lower energies. The fact that some Be/X-ray binaries with major X-ray 
outbursts are located within the positional error boxes of unidentified 
and variable gamma-ray sources (e.g. Zhang et al. 2002; Romero et al. 
2001, 2004) posits the problem of whether such systems are capable 
of generating non-thermal high-energy emission. In this paper we dis- 
cuss how MeV, GeV and even TeV gamma-rays might be produced 
by accreting neutron stars in some Be/X-ray binaries and what is the 
expected time evolution of the high-energy radiation. 
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Spherical accretion 




A transient accretion d\sW is formed 



Figure 1. Sketch of a Be/X-ray binary system 

2. Be/X-ray binaries 

Be/X-ray binaries are the most abundant type of massive X-ray binary 
systems in the Galaxy, with an inferred population in the range 10'^ — 
10^. These systems are formed by a rapidly rotating Be star and a 
neutron star, which is usually observed in an eccentric orbit. In some 
cases, known as Transient X-Ray Sources (TXRS), significant long-term 
X-ray variability is observed. These sources can vary from 1 to 3 orders 
of magnitude going from a quiescent phase to a major outburst, over 
timescales that can go from days to a few years. Sometimes, a major 
outbursts is followed by a series of smaller bursts. Perhaps the most 
extensively studied source of this class is A0535+26 (e.g. Giovanelli &; 
Sabau Graziati 1992), although a few tens are known. 

The general situation in a TXRS is depicted in Figure 1. The Be star 
loses mass in the form of a thick equatorial disk with a low expansion 
velocity and a low-density, fast polar wind. Far from the Be star, the 
neutron star accretes in the Bondi-Hoyle regime. Usually, X-ray flares 
occur once per orbit when the accretion rate increases dramatically 
close to the periastron passage. Occasionally, the equatorial disk of 
the primary can extend beyond the orbit of the compact star, which 
can then penetrate the disk leading to a major outburst. A transient 
accretion disk can be formed around the neutron star in such systems. 
The disk can survive during more than a single orbit in the case of a 
major outburst or it can last only for a fraction of the orbit in the case of 
minor outbursts. Evidence for the existence of such transient disks has 
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Figure 2. Sketch of the magnetosphere of an accreting neutron star when Q,, < fidiak- 
From Romero et al. (2001). 

been found through the detection of quasi-periodic oscillations (QPOs) 
in the power spectra of the hard X-ray flux of some objects (e.g. Finger 
et al. 1996a). 

When the equatorial disk of the Be star retracts a major outburst 
may be followed by a series of minor ones in the same object, as 
observed, for instance, in the system 2S 1417-64 (Finger et al. 1996b). 

3. The magnetosphere of accreting neutron stars 

Once formed, the accretion disk can penetrate the magnetosphere of the 
rotating neutron star (Gosh & Lamb 1979). This penetration creates 
a broad transition region between the unperturbed disk flow and the 
co-rotating magnetospheric flow. In this transition zone, at the inner 
radius Rq the angular velocity is still Keplerian. Between Rq and the 
co-rotation radius -Rco there is a thin boundary layer where the angular 
velocity departs from the Keplerian value. At i?co the disk is disrupted 
by the magnetic pressure and the matter is channeled by the magnetic 
field lines to the neutron star polar cap, where it impacts producing 
hard X-ray emission. 

Cheng & Ruderman (1989, 1991) have studied the magnetosphere 
of accreting neutron stars when the disk rotates faster than the star, 
showing that three definite regions are formed: 1. A region coupled to 
the star by the magnetic field lines that do not penetrate the disk, 
2. An equatorial region linked to the disk by the field attached to it. 
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and 3. A gap entirely empty of plasma separating both regions. In 
this gap E • B 7^ and a strong potential drop Ay is established (see 
Fig. 2 for a sketch of the situation). For typical parameters in Be/X- 
ray binaries, AV ^ 10^^~^^ V (Romero et al. 2001). Protons entering 
into the gap from the stellar co-rotating region will be accelerated to 
multi-TeV energies and directed to the accretion disk by the field lines, 
where they will impact producing tt^ and tt^ in pp interactions with 
the disk material. These pions will decay inside the disk producing 
electron-positron pairs, high-energy (TeV) photons, and neutrinos. The 
neutrinos will escape from the system (see Anchordoqui et al. 2003 
for estimates of the neutrino flux), but the pairs and the gamma-rays 
initiate electromagnetic cascades that will degrade the energy of the 
gamma-ray photons to form a standard cascade spectrum with a cut- 
off at MeV-GcV energies. The Icptons in the disk will lose energy by 
inverse Compton (IC) interactions with the thermal X-ray photons, 
by synchrotron radiation in the field that penetrates the disk, and by 
relativistic Brcmsstrahlung in the ions of the accreting gas. The details 
of the cascade development will be determined by the relative ratios of 
the cooling times among all these processes and the opacities for the 
produced photons in the specific cases. 



4. Gamma-ray production in the accretion disk 

Protons entering into the gap from the stellar co-rotating region are 
accelerated up to energies Ep ~ eVmax, where the potential drop V^ax 
is given by (Cheng Sz Ruderman 1991): 

Vmax - 4 X 10-/?-V2 (^) '''Re'^'Ll^/B^r V. (1) 

Here, ^ 1.4 Mq is the mass of the neutron star, Rq is its ra- 
dius in units of 10^ cm, B12 is the magnetic field in units of 10^^ 
G, and L37 is the X-ray luminosity in units of 10^^ erg s~^. The 
parameter /? = 2ro/rA ~ 1 is twice the ratio of the inner accretion 
disk radius to the Alfven radius (Gosh & Lamb 1979). The Alfven 
radius for spherical accretion can be determined by the condition that 
the unscreened magnetic energy density of the stellar field, B'^/Stt, 
becomes comparable to the kinetic energy density of the accreting 
matter, pv'^/2, i.e. S(rA)^/87r = p{rA_) v{rx)'^/2 , which yields ta = 
/x^/^M~^/^ (2GM*)~^/^, with /j, being the dipolar magnetic moment 
and M the accretion rate. 

The ultra-relativistic proton flux will be directed to the accretion 
disk by the field lines. The maximum current that can fiow through 
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the gap can be determined from the requirement that the azimuthal 
magnetic field induced by tlie current does not exceed that of the initial 
field B: Jmax ~ B^R^r^^ (Cheng k, Ruderman 1989); equivalently 

J^ax ~ 1.5 X 1024/3-2 ~'"rI"l'J,'b^,"' esu s-\ (2) 

The number of protons that impact on the disk per unit of time is 
huge: Np = Jmax/e ~ O (10^^) s"^ and the total power injected in the 
disk results: ^inj ~ JmaxV^ax ~ O (10^^) erg s^^ 

Since all protons arrive to the disk surface with basically the same 
energy, we can estimate the proton flux as: 

' 27r TreRl s cm2 sr ^ ' 

These protons will have a probability of penetrating up to a distance 
L in the disk given by P(L) = 1 - e'^/^^f, with Ap^ = (ap^nf'^^)''^. At 
the relevant energies app ~ 4 X lO^^S cm^. Once the protons interact, 
neutral and charged pions are produced. The decay of the charged pions 
leads to neutrino production in the disk (Cheng et al. 1990, Anchor- 
doqui et al. 2003) and the injection of new e='=-pairs. The gamma-ray 
emissivity from the neutral pions is: 



poo a^lE^) 
«7(£^7) = 2 / . . " , 4 dE,, (4) 



where 



and 



2 4 

ET-{E^) = E^ + ^, (5) 

q^iE,) = J 5{E, - UiE - mp^))app{E)Jp{E) dE 

= T-c^pp (mpC^ + Jp (vripC^ + . (6) 

Jit \ /tt / V /tt / 

Here /jr is the fraction of the kinetic energy of the proton that goes into 
the neutral pion (and then to gamma-rays). Then, using expression (3) 
we obtain: 

q,iE,) ^ llBEi^nAEAsiE, - UEp) . (7) 

e-K j-^Ep atomeVssr 

The intensity of the gamma-ray emission generated in the disk through 
7r°-decays results: 

j-2h 

I°{E^) = q^{E^)E^^npe-'/^^^ dz. (8) 
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Figure 3. Time evolution of the X-ray luminosity of the disk for a major {L^^^^ = 1) 
and a normal outburst (Lg^^*^ = 0.2, dashed curve), t = corresponds to the 
periastron passage. 



The energy released per unit of time in the disk by these photons is 
L^"(ii;^) = 2^i?inj(l-e-2VApp)^ (9) 

and the gamma-ray luminosity emerging from the other side of the disk 
is: 

= L""(^^)e-2'^"^''"''. (10) 

Here, a^p ~ 10^^^ cm^. The gamma-rays that are absorbed can initiate 
electromagnetic cascades inside the disk. 



5. Time-dependent estimates 

In order to calculate the evolution of all relevant parameters during 
both a major and a normal outburst, we have adopted the X-ray light 
curves observed for the Be/X-ray transient 2S 1417-624 (Finger et al. 
1996b), with a normalization^ such that at the peak of the major out- 
burst, which is coincident with the periastron passage, the luminosity 
is L37 = 1. The light curves are shown in Fig. 3. 

In Figs. 4-7 we show the evolution of different parameters that char- 
acterizes the system and the high energy emission for the cases of an 
unscreened disk (r] = 1) and a partially screened disk (rj = 0.2). The 
relation of the screening factor rj and the inner radius of the accretion 
disk is given by (Wang 1996): 

Ro = 1.357o^'r/^/VA, (11) 

^ The distance to this source is not well-established 
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Figure 4- Time-evolution of disk parameters during a major outburst for models 
with 77 = 1 (dashed curve) and 77 = 0.2 (solid curve), t = corresponds to the 
periastron passage. 1. Inner radius of the accretion disk (-Ro). 2. Half-height of the 
accretion disk at _Ro. 3. Magnetic field inside the disk at Rq. 4. Particle density of 
the disk at Rq. 5. Photon density in the disk at Ro- 6. Average energy of the disk 
photons at Rq. 



where 70 = —Bfft^/Bzf, ~ 1 is the magnetic pitch angle. 

In the figures we show, for each case, the evolution of the inner radius 
of the accretion disk, the half-height of the disk, the z-component of the 
magnetic field, the particle density inside the disk, the photon density, 
the average energy of the photons, the energy of the first-generation 
gamma-rays inside the disk, the injected energy into the disk in the form 
of relativistic particles, the opacity of the disk to proton propagation, 
the gamma-ray luminosity produced inside the disk as a consequence 
of the vr'^-decays, the opacity of the disk to these high-energy gamma- 
rays, and the luminosity of the gamma-rays that finally emerge from 
the other side of the disk. All values are given at Rq. 

As we have already mentioned, those gamma-rays that are absorbed 
in the disk can initiate an electromagnetic cascade if the secondary pairs 
are energetic enough as to generate a second generation of gamma-rays. 
In Figure 8 we show the evolution of the different ratios of cooling times 
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Figure 5. Time- evolution of different parameters related to the high-energy emission 
during a major outburst for models with rj — 1 (dashed curve) and rj = 0.2 (solid 
curve), t — corresponds to the periastron passage. 1. Energy of the first-generation 
photons from Tr^-decays. 2. Power injected into the disk in the form of relativistic 
protons. 3. Opacity to proton propagation in the disk. 4. Gamma-ray power from 
TT^-decays produced inside the disk. 5. Opacity to TeV gamma-ray propagation in 
the disk. 6. High-energy gamma-ray luminosity that emerges from the other side of 
the disk. 

for secondary pairs in the disk. During major outbursts in models with 
r/ = 0.2 three generations of pairs are created. The first generation cools 
mainly by synchrotron radiation, whereas in the other generations the 
cooling is dominated by relativistic Bremsstrahlung. The GeV gamma- 
rays produced by the last generation of pairs escape from the disk but 
not necessarily from the system. Inverse Compton losses, which were 
not dominant inside the disk, are now responsible for further cascades 
in the photosphere (see next section). In models of major outbursts 
with r] = 1 only two generations of pairs are produced close to the 
periastron, and just one when the disk is far from it. These pairs cool 
by synchrotron radiation mainly, although Bremsstrahlung takes over 
far in the extremes of the light curve. During minor outbursts cascades 
cannot be sustained: TeV gamma-rays from vr-decays and MeV-GeV 
photons from synchrotron and Bremsstrahlung cooling of pairs origi- 
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Figure 6. Idem Fig. 4, but for a normal outburst. 



nated in the vr -decays are directly injected through the disk into the 
photosphere. 



6. Absorption in the photosphere 

The optical depth for a photon with energy in the photosphere of 
the accreting neutron star is 



/•OO /"OO 

t{Ro,E^)= / n^^{E^,r)a^-^+{E^,E^)drdE^, (12) 

where -Ex is the energy of the ambient thermal X-ray photons, nph(-Ex; f) 
is their density, and ag-^+{Ex, E^) is the photon-photon pair creation 
cross section given by: 



2C(e^-2) + (3-Oln 



1-^ 



(13) 
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Figure 7. Idem Fig. 5, but for a normal outburst. 



where tq is the classical radius of the electron and 



{nieC 



,2^2 



.1/2 



ExE^ 



(14) 



In Eq. 12, -Emin is the threshold energy for pair creation in the X-ray 
field. This field can be considered as formed by two components, one 
from the disk and the other from the neutron star: nph = nph,i + ?T-pii,2- 
Here, 



/7rB{Ex)\ Rl 



with 



B{E^ 



■X, 



he Ex 



2El 



I) 



(15) 



(16) 



(/ic)2 (e-^^x/feTdi3k 

and T^isk = -E'ph.disi^/^) is the disk contribution. The emission from the 
heated polar caps of the star can be approximated by a Bremsstrahlung 
spectrum: 



■nph,2{Ex,r) 



for Ex>1 keV, 



(17) 
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Figure 8. Evolution of the different ratios of cooling times for secondary pairs in an 
average electromagnetic cascade developed in the accretion disk of a systems with 
screening factor rj = 0.2 (dashed) and rj = 1 (solid) during normal (upper panels) 
and major (lower panels) outbursts. 




Log E ^ [eV] Log E ^ [eV] 

Figure 9. Opacity to gamma-ray propagation in the photosphere at the periastron 
passage for both major (solid) and normal (dashed) outbursts and disks with r] = 0.2 
(left) and 77 = 1 (right). 



where Lq is the total luminosity in hard X-rays and kT ^ 100 keV. 

In Figure 9 we show the opacity to gamma-ray propagation at the 
periastron passage as a function of the gamma-ray energy for major 
and normal outbursts in both models with r/ = 0.2 and ry = 1. We 
see that gamma-rays with energies between ~ 1 GeV and ~ 1 TeV 
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Figure 10. Time evolution of the opacity to gamma-ray propagation in the photo- 
sphere for a major outbursts in disks with r) — 0.2 (solid curve) and rj — 1 (dashed 
curve) at two different energies. 
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Figure 11. Idem Fig. fO for a normal outburst. 



energies are absorbed. These gamma-rays initiate inverse Compton 
cascades that degrade their energy. The final spectrum will be soft with 
an index of ~ 2.5 and a high-energy cutoff around ~ 2m1c^ / {Ex) ~ 
250{{Ex) /I keV)~^ MeV (Aharonian et al. 1985, Aharonian & Plyashesh- 
nikov 2003). In Figures 10 and 11 we show the time evolution of the 
opacity at 1 GeV and 1 TeV for different models and outburst types. 
GeV gamma-rays can escape only from minor outbursts and no screen- 
ing of the magnetic field in the disk. Photons with energies above 10 
TeV and below 500 MeV always can escape from the photosphere. 



7. Discussion 



The variability pattern of a given source at gamma rays will depend on 
both the physical parameters that characterize the particular system 
and the energy band of the observations. For instance, almost all GeV 
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sources will be strongly absorbed in their X-ray photosphere during 
major outbursts. In some normal outbursts, however, a flaring GeV 
event might be observable. At TcV energies the gamma rays are par- 
tially absorbed in the disks, especially during major outbursts (a fact 
already noticed by Cheng et al. 1991). The TeV flux can be detectable 
for sources in the inner spiral arms of the Galaxy, once again, mainly 
during normal outbursts. During major flaring X-ray events, the TeV 
gamma rays can be variable or not, depending on the particular combi- 
nation of parameters. At MeV energies the sources should be detectable 
for a wide range of parameters, for both normal and major bursts. A 
correlation between the X-ray and the MeV gamma-ray flux should be 
observed in these systems, whereas an anti-correlation should be the 
case in GcV sources, due to the different absorptions. 

Future detection of Bc/X-ray binaries by GLAST and Cherenkov 
telescopes might be favored by the observation of the candidates during 
normal X-ray outbursts. 
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